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Abstract 

In this paper, we analyze the scalar mesons /q(980) and /q(1500) from the decays B® — > 
/o(980)7T°, /o(1500)7T° within Perturbative QCD approach. From the leading order calculations, 
we find that (a) in the allowed mixing angle ranges, the branching ratio of B® — > /o(980)7r° is 
about (1.0 ~ 1.6) x 10 -7 , which is smaller than that of B® — > fo(980)K° (the difference is a 
few times even one order); (b) the decay B® — > /o(1500)7r° is better to distinguish between the 
lowest lying state or the first excited state for /o(1500), because the branching ratios for two 
scenarios have about one-order difference in most of the mixing angle ranges; and (c) the direct 
CP asymmetries of B® — >■ /o(1500)7r° for two scenarios also exists great difference. In scenario 
II, the variation range of the value A d j P {B° s -> /o(1500)vr°) according to the mixing angle is 
very small, except for the values corresponding to the mixing angles being near 90° or 270°, 
while the variation range of Acp(B® — > /o(1500)7r°) in scenario I is very large. Compared with 
the future data for the decay B® — > /o(1500)7r°, it is ease to determine the nature of the scalar 
meson / (1500). 
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I. INTRODUCTION 



For the underlying structure of the scalar mesons is still under controversy, 
there are two typical schemes for the classification to them [j], 0]. The nonet 
mesons below 1 GeV, including / (600), / (980), if* (800) and a (980), are usually 
viewed as the lowest lying qq states, while the nonet ones near 1.5 GeV, including 
/o(1370),/o(1500)// (1700),A:*(1430) and a (1450), are suggested as the first excited 
states. Here we denote this scheme as scenario I, and the following scheme as scenario II: 
the nonet mesons near 1.5 GeV are treated as qq ground states, while the nonet mesons 
below 1 GeV are exotic states beyond the quark model such as four-quark bound states. 
In order to uncover the inner structures, many approaches are used to research the modes 
of B U) d decaying into a scalar and a pseudoscalar (vector) meson,such as the generalized 
factorization approach [3||, QCD factorization approach (QCDF) |4j-l6|, Perturbative QCD 
(PQCD) approach [7H12||. On the experimental side, along with the running of the Large 
Hadron Collider beauty experiments (LHC-b), some of B® decays involved a scalar in the 
final states might be observed in the Large Hadron Collider beauty experiments (LHC-b) 



13l . [14|. In order to make precision studies of rare decays in the B- meson systems, the 
LHC-b detector is designed to exploit the large number of b hadrons produced. Further- 
more, it can reconstruct a B-decay vertex with very good resolution, which is essential 
for studying the rapidly oscillating B s mesons. Some of B® decays involved a scalar in 
the final states can also serve as an ideal platform to probe the natures of these scalar 
mesons. So the studies of these decay modes for B® are necessary in the next a few years. 

In this paper, we will study the branching ratios and the direct CP asymmetries of 
B® /o(980)7r, /o(1500)7r within Perturbative QCD approach based on kx factorization. 
It is organized as follows: In SectJTTl we introduce the input parameters including the 
decay constants and light-cone distribution amplitudes. In SecJIII[ we then apply PQCD 
approach to calculate analytically the branching ratios and CP asymmetries for our con- 
sidered decays. The final part contains our numerical results and discussions. 



II. INPUT PARAMETERS 

In order to make quantitative predictions, we identify / (980) as a mixture of ss and 
nn = (uu + dd) / a/2, that is 

|/o(980)> = \ss)co$9 + \nn) sinO, (1) 

where the mixing angle 9 is taken in the ranges of 25° < 9 < 40° and 140° < 9 < 165° (l5l |. 
Certainly, /o(1500) can be treated as a qq state in both scenario I and II. We consider that 
the meson / (1500) and /o(980) have the same component structure but with different 
mixing angle. 

For the the neutral scalar meson / (980), / (1500) cannot be produced via the vector 
current, we have (fo(p)\q2l[i<li\ty = 0. Taking the mixing into account, the scalar current 
(/o(p)|?29i|0) = Tn s fs can be written as: 

</o \dd\0) = (fo\uu\0) = -±=m f J%, (f%\ss\0) = m f J s fo , (2) 
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where f^ 1 '^ represent for the quark flavor states for nn and ss components of /o meson, 
respectively. For the scalar decay constants /j? and fj o are very close [5J, we can assume 

ff = fj Q and denote them as ff Q in the following. 

The twist-2 and twist-3 light-cone distribution amplitudes (LCDAs) for different com- 
ponents of fo are defined by: 



(Zo(p)|^)MO),|0> 



1 



V2N< 



c JO 



dxe lxpz {^ fo (x) + m fo $ s fo (x) + m fo {i+i- - l)^ {x)} jh 

(3) 



where we assume fo(p) and fo(p) are same and denote them as fo(p), n + and are 
light-like vectors: n + = (1, 0,0r),n_ = (0,1, Or). The normalization can be related to 
the decay constants: 



dx§f (x 



dx^ T h {x) = 0, 



dx$l (x) = -J^= 

/(,V ' 2y/2W c 



The wave function for it meson is given as [16 

1 



c 



7s [P^{x) + ml^ (x) + Cml^i - v ■ n)^(x)] 



(4) 



(5) 



where P and x are the momentum and the momentum fraction of 7r meson, respectively. 
The parameter ( is either +1 or —1 depending on the assignment of the momentum 
fraction x. 

In general, the B ? meson is treated as heavy-light system and its Lorentz structure can 



be written as [13, [18 



B., 



:^B. + M B MB.{h). 



(6) 



For the contribution of (j>B s is numerically small [19| and has been neglected. 



III. THEORETICAL FRAMEWORK AND PERTURB ATIVE CALCULA- 
TIONS 



Under the two-quark model for the scalar mesons supposition, we would like to use 
PQCD approach to study B° s — > f (980)n, / (1500)7r decays. In this approach, the decay 
amplitude is separated into soft, hard, and harder dynamics characterized by different 
energy scales (t, rriB s , My/)- It is conceptually written as the convolution, 

A(B° S f ir) ~ Jd%d%d 4 k 3 Tr [C(t)^ B .{h)^ f M^(k 3 )H(ki, k 2 , k 3 , t)) , (7) 

where k^s are momenta of anti-quarks included in each mesons, and Tr denotes the 
trace over Dirac and color indices. C(t) is the Wilson coefficient which results from the 
radiative corrections at short distance. The function H(ki, k 2 , k 3 , t) describes the four 
quark operator and the spectator quark connected by a hard gluon whose q 2 is in the 
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order of AMb s , and includes the 0(\JKMb s ) hard dynamics. Therefore, this hard part 
H can be perturbatively calculated. 

Since the b quark is rather heavy, we consider the B° s meson at rest for simplicity. It 
is convenient to use light-cone coordinate {p + ,P~,Pt) to describe the meson's momenta, 

p ± = ^=(p°±p 3 ), and p T = (p 1 ,p 2 ). (8) 

Using these coordinates the B° s meson and the two final state meson momenta can be 
written as 

Pb s = ^|(M,0 T ), P 2 = ^|(1,0,0 T ), P 3 = ^|(0,1,0 T ), (9) 

respectively. The meson masses have been neglected. Putting the anti-quark momenta in 
B®, f and 7r° mesons as k%, k 2 , and fc 3 , respectively, we can choose 

ki = (xxP+^Mt), k 2 = (x 2 P+,0,k 2T ), k 3 = (0,x 3 P 3 -,k 3r ). (10) 

For our considered decay channels, the integration in eq. (JZJ) will lead 

to 

A(Bg —7- fo7t°) ~ J dxidx 2 dx 3 bidbxb 2 db 2 bzdb 3 

•Tr [C(t)$ B . (xi, h)$ fo (x 2 , b 2 )^(x 3 , b 3 )H( Xi , b h t)S t {xi) e" 5(t) ] , (11) 

where b{ is the conjugate space coordinate of k iT , and t is the largest energy scale in func- 
tion H(xi,bi,t). The large logarithms (lnmw/t) coming from QCD radiative corrections 
to four-quark operators are included in the Wilson coefficients C(t). The large double 
logarithms (In 2 Xi) on the longitudinal direction are summed by the threshold resumma- 
tion [20], and they lead to St(xi), which smears the end-point singularities on x^. The 
last term, e~ s ^\ is the Sudakov form factor, which suppresses the soft dynamics effec- 
tively |2l|. Thus it makes the perturbative calculation of the hard part H applicable at 
intermediate scale, i.e., M Bg scale. 

We will calculate analytically the function H(xi,bi,t) for B° s — > f n° decays in the 
leading-order and give the convoluted amplitudes. For our considered decays, the related 
weak effective Hamiltonian H e ff can be written as j22| 



UeH ~ ^ VqbVq * s 



q=u,c 
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(CMOiifj) + c 2 (^)o 2 V)) + E °<M 



i=3 



(12) 



with the Fermi constant Gf = 1.16639 x 10 GeV 2 , and the CKM matrix elements V. 
We specify below the operators in ~H e ff for 6 — > s transition: 



(13) 







up^Lb a , 


U 2 


= Sofl^LUa 


■ up-f^Lbp , 


o 3 


= s a ^Lb a ■ 


Eg' Q.'plM'p > 

Eg' QpinRq'p , 


04 


= Sa^Lbp 


Eg' Qpl^Lq'a , 
Eg' Qpln R Qa , 


o 5 


= s a YLb a ■ 


o e 


= S a ^Lbp 


7 


= \s a ^Lb a 


■ Eg' eq'QpinRqp 

■ Eg' e.q'p^Lq'p 


, Os 


= \s a ^Lb( z 


■ Eg' e g'Q'p7^ R Qa 

■ Eg' e q'Q.'plM'a 


o 9 


= p a YLb a 


010 


= p a ^Lb^ 
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where a and (3 are the SU(3) color indices; L and R are the left- and right-handed 
projection operators with L — (1 — 75), -R = (1 + 75). The sum over q' runs over the quark 
fields that are active at the scale \i = 0(m Bg ), i.e., (q'e{u,d, s,c,b}). 

We will show the whole amplitude for each diagram including wave functions. There 
are 8 type diagrams contributing to the — > f TT° decays are illustrated in FigJTJ We 
first calculate the usual factorizable diagrams (a) and (b). Operators Oi, 2 ,3,4,g,io are (V — 
A)(V — A) currents, and the operators 05,6,7,8 have the structure of (V — A)(V + A), the 
sum of the their amplitudes are written as F e f and F^, respectively. 

F ef = Ff f l = 8TrC F m%Jn / dx x dx 2 / b 1 db 1 b 2 db 2 $B s (xi,bi) 

Jo Jo 

x { [(1 + x 2 )$ fo (x 2 ) - r /0 (l - 2x 2 ) ($ s fo (x 2 ) + $l(x 2 ))] 

xE ei (t)h e (x x , x 2 , 61, b 2 ) - 2r fo <f> s fo (x 2 )E ei (t')h e (x 2 , x u b 2 , b x )} , (14) 

where /„- is the decay constant of 7r meson, r/ = mj /ms s . 

In some other cases, we need to do Fierz transformation for the corresponding operators 
to get right flavor and color structure for factorization to work. We may get (S—P)(S+P) 
operators from (V — A)(V + A) ones. For these (S — P)(S + P) operators, Fig. 1(a) and 
1(b) give 



/»1 /*oo 

Ff f 2 = 167rC F m B J^ / dx x dx 2 \ b x db x b 2 db 2 $ Bs (x 1: bi) 

Jo Jo 

x {- [$ fo {x 2 ) + r fo (x 2 $l(x 2 ) - {x 2 + 2)$%{x 2 ))] 
xE ei (t)/i e (xi,x 2 ,6i,6 2 ) + 2r fo <$> s fo (x 2 )E ei (t')h e (x 2 ,x 1 ,b 2 ,b 1 )} , (15) 

where r n = m^jmB a - 

For the non-factorizable diagrams 1(c) and 1(d), all three meson wave functions are 
involved. The integration of b 2 can be performed using 5 function 5(b 3 — b 2 ), leaving only 
integration of bi and 63. Here we have two kinds of contributions: M e f , and Mjy 2 
describe the contributions from the (V - A) (V - A) , (V - A) (V + A) and (S - P) (S + P) 
operators, respectively. 

f 1 roo 

M efo = 32ixC F m B J \J2N C / dx x dx 2 dx 3 / b 1 db 1 b 3 db 3 ^ Ba (x 1 ,b 1 )^(x 3 ) 

Jo Jo 

x{ - [(x 3 - 1)^/0(^2) -r /o x 2 ($f (x 2 ) - ^ {x 2 ))]E' ei {t)h n {xi,l -x 3 ,x 2 ,h,b 3 ) 
- [(x 2 + x 3 )$ /o (x 2 ) + r /o x 2 ($J o (x 2 ) +&f {x 2 ))] E' ei {t')h n {x l ,x 3 ,x 2 ,b l ,b 3 ) 1,(16) 
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(g) (h) 



FIG. 1: Diagrams contributing to the B® — > /o 71 " decays. 

/*1 /*oo 

Mef = 32irC F m 4 B J \ dx\dx 2 dx 3 / b\db\b 3 db 3 

Jo Jo 

x^xM^ixs - l)$ fo (x 2 )(^(x 3 ) + C(x 3 )) 

+rfo®lM{(x 2 + x 3 - l)$£(z 3 ) + (-x 2 + x 3 - l)$Z{x 3 )) 

+r fo ® S fo (x 2 )((x2 -x 3 + l)$£(z 3 ) - (x 2 + x 3 - l)C(^s))] 
xE' ei (t)h n (x 1 , 1 - x 3 , x 2 , b u b 3 ) + E' ei (t')h n (x 1 ,x 3 , x 2 , 61, b 3 ) 
x[-x 3 $f (x 2 )($Z(x 3 ) - ^(x 3 ))-r fo x 3 ($ s fo (x 2 ) - $J (x 2 )) 

x(^(x 3 )-^(x 3 ))-r f0 x 2 (^%(x 2 ) + ^l(x 2 ))(^(x 3 ) , (17) 
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pi poo 

M p e l = -32nC F m B Jy/2N^ / d Xl dx 2 dx 3 / M&i 63*3 $b.(zi, h)^(x 3 ) 

x| [(^3 - a; 2 - l)$/o(^ 2 ) - r /o x 2 (® S fo (x 2 ) + $/ (x 2 ))] 

xE' ei (t)h n (x 1 , 1 - x 3 , x 2 , 61, 63) + E' ei (t')h n (x 1 ,x 3 , X 2 , 61, 63) 

x [x 2 $ /0 (x 2 ) +r /0 x 2 ($f (x 2 ) - $J (^ 2 ))] } . (18) 

For the non-factorizable annihilation diagrams (e) and (f), again all three wave func- 
tions are involved. For the (V — A)(V — A) and (S — P)(S + P) operators, the results 
are 

pl poo 

M. a f = — 327rCi?m^ s / a/ 2N c / dx±dx 2 dx 3 / bidbib 3 db 3 

Jo Jo 

^ Bs { x iM)\yE' ai (t)h na (x 1 ,x^x 2 MM) [x3&*(xa)&f (x 2 ) 

+rnr fo ^ (x 2 )((x 2 -x 3 + l)$l(x 3 ) - (x 2 + x 3 - l)Q*(x 3 )) 
+r«r h <$>%{x 2 ){{-x 2 + x 3 + 3)^(x 3 ) + (x 2 + x 3 - l)$?(x 3 ))] 

+Ki(t')h'na(^1^3,X 2 ,b 1 ,b 3 ) [{X 2 - l)^(x 3 )$ /o (x 2 ) 

+r 7r r fo ^ o (x 2 ) ((-x 2 +x 3 + l)$%(x 3 ) - (x 2 + x 3 - \)^{x 3 )) 
+ Vfo <£> s fo (x 2 ) ((x 2 -x 3 -l)^(x 3 ) + (x 2 + x 3 -l)^(x 3 ))] ) , (19) 



/*1 POO 

M^l = -32nC F m%/y/2N^ / dx x dx 2 dx 3 \ M&i M&3 $b.(si, &i) 

Jo Jo 

x I [(x 2 - l)<f> fo (x 2 )<f>*(x 3 ) - 4r n r fo <S> s fo (x 2 )<S>*(x 3 ) + r n r fo ((x 3 - x 2 - 1) 

x (^(a*)** (x 2 ) + <S>T(x 3 )$l(x 2 )) - (x 2 +x 3 -l) {^(x 3 )^ o (x 2 ) 
-$l(x 3 )$ s fo {x 2 )))] E' ai (t)h na ( Xl ,x 3 , x 2 , b u b 3 ) + E' ai (t')h' na ( Xl ,x 3 , x 2 , 61, 63) 
x [x 3 <£> fo (x 2 )^(x 3 ) + x 3 r n r fo ($ s fo (x2) - $J(x2))(^(x s ) + ^(x 3 )) 

+r w r ib (l-x 2 )($J ) (x2)+^(x2))(^(xs)-^(xs))] } . (20) 

The factorizable annihilation diagrams (g) and (h) involve only the n and fo mesons' 
wave functions. There are three kinds of decay amplitudes for these two diagrams. F a f 
is for (V — A)(V — A) type operators, F^} is for (V — A)(V + A) type operators, while 
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F^ o is for (S — P)(S + P) type operators: 

/l POO f 

dx 2 dx 3 J b 2 db 2 b 3 db 3 | [(x 2 - l)$^(x 3 )$ /o (x 2 ) 

+2r n r fo (x 2 - 2)^(x 3 )$%(x 2 ) - 2r„r fo x 2 <S>? {x 3 )$ T fo (x 2 )] 

xE ai (t)h a (x 3 , 1 - x 2 , b 3 , b 2 ) + E ai (t')h a (l - x 2 , x 3 , b 2 , b 3 ) [x 3 $^(x 3 )$ fo (x 2 ) 

+2r n r f0 ^%(x 2 )((x 3 + l)^(x 3 ) + (x 3 - 1)*J(* 8 ))] 1, (21) 



/»1 /*oo 

= 16irC F m Ba f Ba / cfe 2 <ir 3 / b 2 db 2 b 3 db 3 
Jo Jo 

x {[r/ (x 2 - l)<^(x 3 )($* (x 2 ) + $J o (x 2 )) - 2r^(x 3 )$ fo (x 2 )} 
xE ai (t)h a (x 3 , 1 - x 2 ,b 2 ,b 3 ) - E ai {t')h a {\ - x 2 ,x 3 ,b 2 ,b 3 ) 
x[2r /0 $^(x 3 )<l>f (x 2 ) + r A $/ (a ;2 )(^N - (22) 

If we exchange the 7r° and /o in FigJT], the result will be different. In the considered 
decays, the meson f cannot lie in the emitted position, like 7r' position in Fig.l(a-d). So 
only the annihilation type diagrams left, just like Fig.l(e-h), can give contributions. They 
are listed as follows: 

Man = 32nC F m B J^/2N c J dx x dx 2 dx 3 J M&i b 2 db 2 $b.{xi, h)l [-x 2 $£(x 3 ) 

x$ /o (x 2 ) + r n r fo $ T fo (x 2 ) ({x 2 + x 3 - l)^(x 3 ) + {-x 2 + x 3 + l)$£(ar 3 )) 
+r w r fo <S> s fo (x 2 ) ((x 2 -x 3 + 3)^(x 3 ) - (x 2 + x 3 - l)$£(x 3 ))] 
xE' ai (t)h na (x 1 , x 2 , x 3 , b u b 2 ) - E' ai (t')h' na (xi, x 2 , x 3 , b u b 2 ) [(x 3 - l)$^(x 3 ) 
x$ /o (x 2 ) + r v r h $%(x 2 ) {{x 2 - x 3 + l)$£(x 3 ) - (x 2 + x 3 - l)$£(z 3 )) 

+r n r f ^l(x 2 ) {(x 2 + x 3 - 1)$; (x 3 ) - (1 + x 2 - x 3 )^(x 2 ))] ) , (23) 



M p a l = -32nC F m B J^/2N^ dx x dx 2 dx 3 b x db x b 2 db 2 ® Bs { x u h)i [±V h $%(x 2 ) 

x$£(x 3 ) + (x 3 - l)$ /0 (x 2 ) $£(x 3 ) + r % r fo {{x 2 -x 3 -l) (x 3 )$ 6 fo (x 2 ) 
-<S>l(x 3 )$ T fo (x 2 )) -{x 2 + x 3 - l)(^(x 3 )$ T fo (x 2 ) - C(x 3 )$J (x 2 ))] 
xE'^Ka^x^XiMM) + [x2®f (x 2 )®i(x 3 ) - x 2 r^r fo (^ s fQ {x2) + $^(x2)) 
x(^(x 3 ) - $T(x 3 )) - r„r fo {l - x 3 ){$>%{x2) - $J o (x2))(^(x 3 ) + ^(x 3 ))] 

xE' a ^t')h' na {x u x 2 ,x 3 MM)\ , (24) 
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/I /"OO ( 

dx 2 dx 3 J b 2 db 2 b 3 db 3 | [(x 3 - l)$^(x 3 )$ /o (x 2 ) 

-2r.r /o (x 3 - 2)<^(x 3 )$f (x 2 ) + 2r 7r r /o x 3 C(^)$/ (x 2 )] 
xE ai (t)h a (x 2 , 1 - x 3 , 6 2 , h) + E ai (t')h a (l - x 3 , x 2 , h, b 2 ) 



x 



[x 2 ^(x 3 )$ /0 (x 2 ) - 2r 7r r /o <&£(x 3 )((x 2 + 1)<&J (x 2 ) + (x 2 - l)Sj,(s 2 ))] , (25) 



i>1 poo 

Fan = -^C F m B J Bs / dx 2 dx 3 / b 2 db 2 b 3 db 3 

Jo Jo 



x j[rv(x 3 - l)$ /o (x 2 )(^(x 3 ) +^;(x 3 )) + 2r fo $ 7r (x 3 )&} o (x 2 )] 
xE ai (t)h a (x 2 , 1 - x 3 ,b 2 ,b 3 ) - E ai (t')h a (l - x 3 ,x 2 ,6 3 ,6 2 ) 
xpr^^^^+r^x^^xg)^^)-^^))]! . (26) 

In the above formulae, the function E are defined as: 

E ei (t) = a s (t) exp[-S B (t) - S 3 (t)], (27) 

E' ei (t) = a s (t) exp[-5 fl (t) - S 2 (t) - S 3 (t)]\ b3=bl , (28) 

£«i(t) = a s (t) exp[-S 2 {t) - S 3 (t)}, (29) 

KM = <x*(t) exp[-S B (t) - S 2 {t) - S 3 {t))\ b3=b2 , (30) 

where a s is the strong coupling constant, S is the Sudakov form factor. In our numerical 
analysis, we use the one-loop expression for the strong coupling constant; we use c = 0.4 
for^the parameter in the jet function. The explicit form of h and S have been given in 

Combining the contributions from different diagrams, the total decay amplitudes for 
these decays can be written as: 

M(B° S -»■ f ir) = M sS x cosfl + —M n nSm9, (31) 

v2 

where 9 is mixing angle, and 

M s -s = V ub V: s (F ef0 a 2 + M ef0 C 2 ) - h tb V t * s [F e p *(a 9 - a?) + M%C W + M^C 8 ] , (32) 

M nn = V ub V: s [(M w + M af0 )C 2 + (F a7T + F af0 )a 2 ] - ^V tb V* [{M m + M afo )C 10 

+{M^ + M2)C* + (F a7r + F af0 )(a 9 - or)] , (33) 



where the combinations of the Wilson coefficients are defined as usual 24, 25 



a% = C 2 + d/3, a 3 = C 3 + C 4 /3, a 5 = C 5 + C 6 /3, a 7 = C 7 + C 8 /3, a 9 = C 9 + C 10 /3, (34) 
a 2 = d + C 2 /3, a 4 = C 4 + CaA a 6 = C 6 + C B /3, a 8 = C 8 + C 7 /3, a 10 = C 10 + C B /3. (35) 
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IV. NUMERICAL RESULTS AND DISCUSSIONS 



The twist-2 LCDA $/ can be expanded in the Gegenbauer polynomials: 



j oo 

$ /0 (x,/i) = -^^/ /o ( /U )6x(l-x)^5 m (/i)C^ 2 (2x-l), (36) 
^V^v c m=l 

where B m (/j) and Cm 2 (a;) are the Gegenbauer moments and Gegenbauer polynomials, 
respectively. The values for Gegenbauer moments and the decay constants are taken (at 
scale \x = lGeV) as: 

Scenario I : / /o(980) = (0.37 ± 0.02)GeV, / /o(1500 ) = -(0.255 ± 0.03)GeV, 

Bi(980) = -0.78 ±0.08, 5 3 (980) = 0.02 ± 0.07, 

Bi(1500) = 0.80 ± 0.40, 5 3 (1500) = -1.32 ± 0.14; 

Scenario II : / /o(150 o) = (0.49 ± 0.05)GeV, 

^(1500) = -0.48 ±0.11, B 3 (1500) = -0.37 ±0.20. (37) 

As for the explicit form of the Gegenbauer moments for the twist-3 distribution am- 
plitudes $j o and $J o , although they have been studied in the Ref. |26|, we adopt the 
asymptotic form: 

^ = 27W/*- * T >» = 27mM l - 2l) - (38) 

The twist-2 pion distribution amplitude and the twist-3 ones <3>f an d ®n nave been 
parametrized as 

-0.028(1 - 14(1 - 2xf + 21(1 - 2x) 4 )] , (39) 

-0.11/8(3 - 30(1 - 2a;) 2 + 35(1 - 2x) 4 )] , (40) 
$Z(x) = ^=(1 - 2x) [1 + 0.56(1 - lOx + 10x 2 )] . (41) 

The B s meson's wave function can be written as: 

M 2 x 2 1 

<P Ba (x,b) = N Ba x 2 (l -x) 2 exp[-^- - ~{u bs b)\ (42) 

where Ub B is a free parameter and we take Ub s = 0.5 ±0.05 GeV in numerical calculations, 
and Nb s = 63.67 is the normalization factor for Ub a = 0.5. 

For the numerical calculation, we list the other input parameters in Table I. 
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TABLE I: Input parameters used in the numerical calculation 



27] 



Masses 


m Mm = 0.98 GeV, m£ = 1.3 GeV, 
"i/o(i500) = 1-5 GeV M Bs = 5.37 GeV, 


Decay constants 


f Bs = 0.23 GeV, A = 0.13 GeV, 


Lifetimes 


t b o = 1.466 x 10~ 12 s, 


CKM 


V tb = 1.0, V ts = 0.387, 
V us = 0.2255, V ub = 0.00393e- i60 ° . 



TABLE II: Decay amplitudes for 5° ->■ / (980)7r°, / (1500)7r° (xlO" 2 GeV 3 ). 



ss 


Fl> T 


K17T 

± e 


mV 1 ' 


Ml 


/ (980)vr° (SI) 


8.98 


1.50 


-3.51 + 5.75i 


-0.01 + 0.04z 


/o(1500)7T° (SI) 


-11.9 


-1.03 


-5.89 + 4.70f 


-0.02 + 0.04i 


/o(1500)^° (SII) 


17.8 


2.38 


-0.77 + 3.64i 


0.001 + 0.02z 


nn 


Ml°> T + MS' T 


M£ + Ml 


Fi°' T + FZ' T 


Fl° + Fl 


/ (980)(nn)vr° (SI) 


6.3 + 6.9i 


-0.05 + 0.03? 


0.89 - O.lli 


-0.001 + 0.025i 


/o(1500)(nn)7r° (SI) 


-13.8 + 20.8i 


0.11 + 0.03i 


-0.015 + 1.34? 


-0.02 - 0.017i 


/o(1500)(nn)^° (SII) 


15.3 -3.1i 


0.01 - O.Oli 


1.62 - 0.63i 


0.04 + 0.04i 



If /o(980) and /o(1500) are purely composed of nn(ss), the branching ratios of B® — > 
/o(980)7r°,/ (1500)7r are: 

B{B° S / (980)(nn)7T°) = (0.461™;™) x 10" 8 , Scenario I, (43) 

B(B° S / (1500)(nn)7T°) = (2.41°;^^) X 10" 8 , Scenario I, (44) 

B{B° S ->• / (1500)(nn)7r°) = (l.l±C;f±^;|±g;f ) x lO" 8 , Scenario II; (45) 

B(B° S /o(980)(ss)vr°) = (1.7^1^1^) x 10~ 7 , Scenario I, (46) 

B(B° S / (1500)( S s)7T°) = (O.eelHo-ais-aos) x lO" 7 , Scenario I, (47) 

B{B° S / (1500)(ss)tt ) = (3.8^1°;^) x 10" 7 , Scenario II, (48) 

where the uncertainties are from the decay constant of fo, the Gegenbauer moments Bi 
and -B3. One can see that the values of B{B S — > fo{nn)7T°) are smaller than the correspond- 
ing values of B{B S / (ss)7r°), it is contrary to the case of B s -> f (980)K°, f (1500)K° 
decays (28| . 

In Table II, we list values of the factorizable and non-factorizable amplitudes from the 
emission and annihilation topology diagrams of B° s — > / (980)7r°, / (1500)7r°. a \ and 

Af^A are the 7r° emission (annihilation) factorizable contributions and non-factorizable 
contributions from penguin operators respectively. Similarly, F ? , and M ? , denote as 
the contributions from fo emission (annihilation) factorizable contributions and non- 
factorizable contributions from penguin operators respectively. F^? 1, , Mj^y denote 
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FIG. 2: The dependence of the branching ratios for 5° -> /o(980)vr° (a) and 5° -> / (1500)7r° 
(b) on the mixing angle using the inputs derived from QCD sum rules. The vertical bands 
show two possible ranges of 0: 25° < 9 < 40° and 140° < 6 < 165°. For the right panel, the 
solid (dotted) curve is plotted in scenario II (I). 



the corresponding contributions from tree operators 02,0\. From the table, regard- 
less of the CKM suppression, one can find that the contributions from tree operators are 
(much) larger than the corresponding ones from penguin operators, especially for the non- 
factorizable emission diagrams and the annihilation diagrams. In fact, the tree operators 
contributions are strongly CKM-suppressed compared to penguin operators contributions. 
Here considering the amplitudes of annihilation diagrams is necessary. 

In Fig. [2J we plot the branching ratios as functions of the mixing angle 9. Because there 
are many discrepancies about the mixing of quark components for the meson / (1500), 
we show the dependence of the branching ratio for — > / (1500)7r on all the mixing 
angle values, that is (0°,360°). In the allowed mixing angle ranges, the branching ratio 
of B° s /o(98GV is: 

1.0 x KT 7 < B(B° S -»■ /o(980)7r°) < 1.6 x 1(T 7 , (49) 

which is smaller than the branching ratio of B® — > fo(980)K°. The difference is a few 
times even one order. As to the decay B® — » / (1500)7r°, it is interesting that this channel 
is better to distinguish between the first excited state (scenario I) and the lowest lying 
state (scenario II) for / (1500). The branching ratio of B° s — > / (1500)7r° for scenario I is 
at the order of 10~ 8 in (0°, 360°), while its value for scenario II is at the order of 10~ 7 in 
most of mixing angle ranges, except for (60°, 120°) and (240°, 300°). If the mixing angle is 
not in these two ranges, it is ease to determine the nature of /q(1500). If the observation 
of the decay B° s — > f (1500)ir at the order of 10~ 7 , it would indicate that scenario II 
is favored. We also find the branching ratios of these decays are smaller than those of 
corresponding channels B° s — > f (980)K°, f (1500)K° about a few times even one order. 

Now we turn to the evaluations of the CP-violating asymmetries of the considered 
decays in the PQCD approach. For the neutral decays B° s — > / (980)7r°, / (1500)7r°, there 
are both direct CP asymmetry A^j r p and mixing-induced CP asymmetry A™p . The time 
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TABLE III: Direct CP asymmetries (in units of %) of -> / (980)vr , / (1500)vr° decays for 
nn and ss components, respectively. 



Channel 


Scenario I 


Scenario II 


Bs -> /o(980)(nn)vr 


51.5 




B° s -> /o(980)( S s)7r° 


-16.0 




£ s ° -> /o(1500)(nn)7r° 


29.0 


15.1 


B° s /o(1500)( S s)7r° 


19.3 


-6.7 




FIG. 3: The dependence of the direct CP asymmetries for B® — > /o(980)7r° (a) and B® — »• 
/o(1500)7T° (b) on the mixing angle 6. For the right panel, the solid (dotted) curve is plotted in 
scenario II (I). The vertical bands show two possible ranges of 8: 25° < 8 < 40° and 140° < 9 < 
165°. 



dependent CP asymmetry of B s decay into a CP eigenstate / is denned as: 

A CP (t) = Af, r P (B s -> /) cos(Am s ) + A^(B S -> /) sin(Am s ), (50) 

with 

I \ 12 _ i o j \ 

Af P {B s /) = ^( 5s /) = _____ (51) 

where 77 = ±1 depends on the CP eigenvalue of /, Am s is the mass difference of the two 
neutral B s meson eigenstates. Here we only give the direct CP-violating asymmetry. 

The direct CP asymmetries of B® — > f (nn)7r°, fo(ss)7r° are listed in Table III. From 
the definition of the direct CP asymmetry Eq. fl5ip and Eq. (j52j) . we can find the sign of 
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Afj r P is determined by the formula: 



(F:> T + M?f) - Mlf {Fl + M- r ) , for/o (aa) ; 
Mj Mf - MjM? = \ {Mlf + Mff + Fff + F£?) (m/° + + F* + f£ ) - (m# 

+Kf + <f + <f ) {Mt + M: >r + Fj° + Fl r ) , for/ (nn) . 



If M^Mf — Mf > 0, the sign of the corresponding direct CP asymmetry is positive, 
contrarily, the value of Afj r P is minus. So one can understand that though the penguin 
operators contributions are much smaller than the tree operators contributions (the differ- 
ence is about two or three orders, seen in Table II), they are important to determine the 
direct CP asymmetry From Table III, it is found that A% r P (B° / (1500)(»7r ) > in 
scenario I, which is contrary to the sign in scenario II. The direct reason is that there ex- 
ists an opposite sign for the 7r° emission factorizable contributions from penguin operators 
between these two scenarios. We also find the direct CP asymmetries for nn components 
of / (980) and / (1500) are larger than those for ss components. 

In Fig. 3, we plot the direct CP asymmetries of the decays B° s — > / (980)7r° and — > 
/o(1500)7T°. One can see the direct CP asymmetry of B° s — > / (980)7r° is: 

- 30% < A d j P (B° s -> / (980)7r°) < -25%, (53) 

for the mixing angle 9 in the range of 25° < 9 < 40°. While if the mixing angle 9 is taken 

in the range (140°, 165°), the value of A% r P {B° s ->• / (980)vr ) is about (-12 5)%. For 

the decay B Q S — > /o(1500)7r°, if the parameters in scenario II are used, one can find the 
variation range of Afj r P {B^ — > /o(1500)7r°) according to most of the mixing angles is very 
small, except for the values for mixing angles near 90° or 270°. while in scenario I, the 
variation range of A^p{B° s — > / (1500)7r°) is very large. The great differences of decay 
constant and Gegenbauer moments of /o(1500) result that there exists great difference for 
the direct CP asymmetries in two scenarios. It gives the hint that one can determine the 
nature of the meson / (1500) by comparing with the future experimental values for these 
direct CP-violating asymmetries. 

It is noticed we consider that the meson /o(1500) is dominated by the quarkonium 
content, the detail discussion can be found in [28]. If we take the mixing mechanism for 
/o(1500) as | / (1500)} = -0.54 | nn) + 0.84 | ss) + 0.03 | G) H and neglect the small 
component of glueball, we have: 

B(B° S /o(1500)vr°) = (4.46l°;^f 1°;!) x 10" 8 , Scenario I, 
B(B° S /o(1500)7T°) = (2.81^;™^) x 10 - 7j Scenario n . 

A d S r P (B° s / (1500)7r°) = -(27.5i°;H:^: 3 )%, Scenario I, 

A d S r P {B° s /o(1500)7T°) = -{9.7t° ° t\itli)%, Scenario II, (54) 

which are the values corresponding to the mixing angle 327.3°. The uncertainties are from 
the decay constant of fo, the Gegenbauer moments B\ and B3 for twist-2 LCDAs of the 
scalar mesons. Certainly, it is only the leading order results. In this process, the 7r° emis- 
sion factorizable amplitudes from the tree operators correspond to the color-suppressed 
tree amplitudes, which are known to be modified by the inclusion of the next-to-leading- 
order (NLO) corrections. From the calculations of the partial NLO corrections j25j, our 



14 



argument is that the NLO contributions might have a small influence on the branching 
ratio. But it is difficult to say that the predicted discrepancy in the CP asymmetries must 
hold under all of the NLO corrections. 

V. CONCLUSION 

In this paper, we study — > f (980)ir°, / (1500)7r° decays in the PQCD factorization 
approach and calculate their branching ratios and the direct CP-violating asymmetries. 
Several remarks are in order: 

• If /q(980) and /o(1500) are purely composed of nn or ss, one can see that the values 
of B(B S — > f (nn)ii ) are smaller than those of B(B S — > f (ss)ii ), it is contrary to 
the B s -> / (980)ir°, / (1500)K° decays. 

• In the allowed mixing angle range, the branching ratio of B° s — > / (980)7r° is: 



which is smaller than that of the decay B° s — > f (980)K . The difference is a few 
times even one order. 

• The decay B° s — > / (1500)7r° is better to distinguish between the lowest lying state or 
the first excited state for / (1500). Because its branching ratios for the two scenarios 
have about one order difference in most of the mixing angle ranges. For example, if 
we take the mixing mechanism for /o(1500) as | /o(1500)} = —0.54 | nn) + 0.84 | ss), 
which corresponds to the mixing angle taking about 327.3°, one can find 



• There also exists great difference for the direct CP asymmetries of the decay B° s — > 
/o(1500)7T° in two scenarios. If the parameters in scenario II are used, one can 
find the variation range of the value Aq P (B^ — > / (1500)7r°) according to most 
of the mixing angles is very small, except for the values corresponding to mixing 
angles being near 90° or 270°, while in scenario I, the variation range of Aqp(B® — > 
/o(1500)7T°) is very large. Certainly, the NLO contributions may give these direct 
CP asymmetries some corrections. 

Acknowledgment 

This work is partly supported by Foundation of Henan University of Technology under 
Grant No. 150374. The author would like to thank Hai-Yang Cheng, Cai-Dian LU, Wei 
Wang, Yu-Ming Wang for helpful discussions. 



[1] N.A. Tornqvist, Phys. Rev. Lett. 49, 624 (1982). 



1.0 x HT 7 < B(B° S -> /o(980)7r°) < 1.6 x 10~ 7 , 



(55) 



B(B° S /o(1500)tt ) = 4.46 x 10~ 8 , Scenario I, 
B(B° S /o(1500)tt ) = 2.81 x 10" 7 , Scenario II. 



(56) 



15 



[2] G.L. Jaffe Phys. Rev. D 15, 267 (1977); Erratum-ibid.Phys. Rev. D 15 281 (1977); 
A.L. Kataev, Phys. Atom. Nucl. 68, 567 (2005), Yad. Fiz. 68, 597(2005); A. Vijande, 
A. Valcarce, F. Fernandez and B. Silvestre-Brac, Phys. Rev. D72, 034025 (2005). 
A. K. Giri , B. Mawlong, R. Mohanta Phys. Rev. D 74, 114001 (2006). 
H. Y. Cheng, K. C. Yang Phys. Rev. D 71, 054020 (2005). 
H. Y. Cheng, K. C. Yang, Phys. Rev. D 73, 014017 (2006). 
H. Y. Cheng, C.K. Chua, K. C. Yang, Phys. Rev. D 77, 014034 (2008). 
C. H. Chen, Phys. Rev. D 67, 014012 (2003). 
C. H. Chen, Phys. Rev. D 67, 094011 (2003). 

W. Wang, Y. L. Shen, Y. Li, C. D. Lii Phys. Rev. D 74, 114010 (2006). 
Y. L. Shen, W. Wang, J. Zhu and C. D. Lii, Eur.Phys.J.C 50:877-887 (2007). 
Z.Q. Zhang and Z.J. Xiao, Chin.Phys.C33(07):508-515 (2009). 
Z.Q. Zhang and Z.J. Xiao, Chin.Phys.C34(05):528-534 (2010). 

N.Brambilla et al, (Quarkonium Working Group), CERN-2005-005, |hep-ph/0412158[ M.P. 
Altarelli and F.Teubert, Int. J. Mod.Phys. A 23, 5117 (2008). 

M. Artuso et al, "B, D and K decays", Report of Working Group 2 of the CERN workshop 
on Flavor in the Era of the LHC, Eur. Phys. J. C 57:309-492 (2008). 
HY. Cheng, Phys. Rev. D 67, 034024 (2004). 
P. Ball and R. Zwicky, Phys. Rev. D 67, 014015 (2005). 

A.G.Grozin and M.Neubert, Phys. Rev. D 55, 272 (1977); M.Beneke and T.Feldmann, 
Nucl. Phys.B 592, 3 (2001). 

H.Kawamura, J.Kodaira, C.F.Qiao and K.Tanaka, Phys. Lett. B 523, 111 (2001); 
Mod.Phys.Lett.A18, 799 (2003). 
C.D.Lii, M.Z.Yang, Eur.Phys.J.C 28, 515 (2003). 
H. N. Li, Phys. Rev. D 66, 094010 (2002). 
H. N. Li and B. Tseng, Phys. Rev. D 57, 443 (1998). 

G. Buchalla, A. J. Buras, M. E. Lautenbacher, Rev. Mod. Phys. 68, 1125 (1996). 
Y. Y. Keum, H. N. Li, A. I. Sanda, Phys. Lett. B 504, 6 (2001); Phys. Rev. D 63, 054008 

(2001) ; C. D. Lii, K. Ukai, M. Z. Yang, Phys. Rev. D 63, 074009 (2001); Y. Y. Keum, H. 
N. Li, Phys. Rev. D 63, 074006 (2001); C. D. Lii , M. Z. Yang, Eur. Phys. J. C 23, 275 

(2002) . 

[24] A. Ali, G. Kramer, CD. Lii, Phys. Rev. D 58, 094009 (1998). 
[25] Z.J. Xiao, Z.Q. Zhang, X. Liu, L.B. Guo, Phys. Rev. D 78, 114001 (2008). 
[26] CD. Lii, Y. M. Wang and H. Zou, Phys. Rev. D 75, 056001 (2007). 
[27] Particle Data Group, C. Amsler et al, Phys. Lett. B 667,1 (2008). 
[28] Z.Q. Zhang, J. Phys. G 37, 085012 (2010). 

[29] H.Y. Cheng, C.K. Chuan and K.F. Liu, Phys. Rev. D 74, 094005 (2006). 



16 



